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1 "Borehole Surveying" 

2 

3 This invention relates to a method and apparatus for 

4 use in surveying of boreholes. 
5 

6 It is known in directional drilling, for example, to 

7 detect the orientation of a drillstring adjacent to 

8 the bit by means of a sensor package for determining 

9 the local gravitational [GX,GY,GZ] and magnetic 

10 [BX,BY,BZ] field components along mutually 

11 orthogonal axes, and to derive from these the local 

12 azimuth (AZ) and inclination (INC) of the 

13 drillstring. Conventionally, the measurements are 

14 made by providing within the instrument package 

15 three mutually perpendicular accelerometers and 

16 three mutually perpendicular magnetic fluxgates. 
17 

18 The present invention is concerned with an 

19 arrangement which requires only two measurement 

2 0 devices, namely a single accelerometer and a single 

21 magnetic fluxgate or a single accelerometer and a 

22 single rate gyro, the latter being preferred for 

23 situations in which magnetic interference is likely 

24 to be encountered. 

25 ' - 

2 6 Accordingly, the present invention provides a method 

2 7 of surveying boreholes, comprising: 



2 

1 providing an instrument package in the leading 

2 end of a drillstring, the instrument package 

3 comprising first and second single-axis sensors 

4 mounted for rotation with the drillstring about the 

5 rotational axis of the drillstring, the first sensor 

6 being an accelerometer and the second sensor being a 

7 magnetic fluxgate or a rate gyro; 

8 v rotating the drillstring; 

9 deriving from the first sensor the inclination 

10 angle of the drillstring at the instrument package ; 

11 and 

12 deriving from the second sensor the azimuth 

13 angle of the drillstring at the instrument package. 
14 

15 Each of the sensors will typically be positioned in 

16 one of. two configurations. In the. first 

17 configuration, the sensor is radially spaced from 

18 the borehole axis and has its sensing axis in a 

19 plane containing the borehole axis and an axis 

2 0 perpendicular thereto. In the second configuration, 

21 the sensor is radially spaced from the borehole axis 

22 and has its sensing axis in a plane parallel with 

23 the borehole axis. 
24 

25 Preferably, the drilling control rotation angle is 

26 also obtained from the sensor outputs. 

27 ... 

2 8 Preferably, the sensor outputs are integrated over 

2 9 the four quadrants of rotation and the desired 

3 0 output angle is derived from the integrated output. 

31 The instrument package suitably includes rotation 

32 angle reference means for use in the integration. 



3 

1 * 

2 Additional information may be derived, such as the 

3 local gravitational and magnetic field" vectors. 
4 

5 . From another aspect, the invention provides 

6 apparatus for use in surveying boreholes, the 

7 apparatus comprising an instrument package adapted 

8 to be included in the leading end of a drillstring, 

9 the instrument package comprising first and second 

10 single-axis sensors mounted for rotation with the 

11 drillstring about the rotational axis of the 

12 drillstring, the first sensor being an accelerometer 

13 and the second sensor being a magnetic fluxgate or a 

14 " rate gyro; and computing means for deriving from the 

15 first sensor while the drillstring is rotating the 

16 inclination angle of the drillstring at the 

17 instrument package, and for deriving from the second 

18 sensor while the drillstring is rotating the azimuth 

19 angle of the drillstring at the instrument package. 
20 

21 The computing means preferably operates to integrate 

22 the sensor outputs over the four quadrants of 

23 rotation and to derive the desired output angle 

24 from the integrated output. 
25 

2 6 The apparatus may further include rotation angle 

2 7 reference means for use in the integration. 

28 . 

2 9 Examples of the present invention will now be 

30 described, by way of illustration only, with - 

31 reference to the drawings, in which: 
32 



4 

1 Fig. 1. illustrates, in general terms, the 

2 operation of a single axis sensor in a drillstring 

3 ■ for sensing any given vector V; 

4 Fig. 2 is a block diagram of one circuit which 
5' may be used to identify rotation quadrant;" 

6 Fig. 3 illustrates the operation where the 

7 sensor is an accelerometer ; 

8 Fig. 4 illustrates the operation where the 

9 sensor is a fluxgate; 

10 Fig. 5 illustrates the derivation of azimuth 

11 angle; and 

12 Fig. 6 illustrates the operation where the 

13 sensor is a rate gyro. 
14 

15 

16 Single -axis sensor 

18 The operation of a single-axis sensor in a drill 

19 string will first be described in general terms. 

20 The application of this to specific sensors is 

21 discussed below. 
22 

23 Referring to Fig. 1, a single-axis sensor 10 is 

24 mounted on a drill string (not shown) . The sensor 

25 10 senses a fixed vector {v} and is mounted in one 

26 of two configurations. 
27 

28 In the first configuration, the sensor 10 lies in a 

2 9 plane containing the rotation axis (OZ) of the drill 

3 0 string and axis (OX) perpendicular to (OZ) . Axis 

31 (OY) makes up the conventional orthogonal set of 

32 axes [OX,OY,OZ] . ; . The sensor 10 is mounted at a 



5 

1 distance r from the (OZ) axis and the angle between 

2 the sensing axis (OS) and the rotational axis (OZ) 

3 is m. 
4 

5 In the second configuration, the sensor 10 is 

6 . mounted in a plane which is parallel to the borehole 

7 axis (OZ) and with its sensing axis perpendicular to 

8 the axis (OY) and making angle m with the direction 

9 of the borehole axis (OZ) . 
10 

11 If the rate of rotation about the (OZ) axis is w and 

12 the components of {v} are {VOZ} along the (OZ) axis 

13 direction and {VOXY> in the (OXY) plane, then if the 

14 output from the sensor 10 for both configuration 1 

15 and configuration 2 of Figure 1 is of the form 
16 

17 V(t) = VOZ.cos(m) + VOXY.sin(m) .cos(w.t) + c 

18 

19 where time t = 0 when the axis (OX) is coincident 

2 0 with the direction of {VOXY} and c is constant for 

21 any fixed rotation rate w. 
22 

23 Thus, the sensor output at time t can be written: 
24 

25 V(t) = Kl.cos(w.t) + K2 . . (i) 

26 

2 7 where Kl = VOXY. sin (m) and K2 = VOZ. cos (m) + c are 

2 8 constant if the vector amplitudes VOZ and VOXY are 

2 9 constant . 
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1 Sensor output integration 
2 

3 The integration of V(t) from any initial time ti to 

4. ti + T/4, where T = 2.7t/w, the time for one 

5, revolution about (OZ) , is 
6 

ti+T/4 ti+T/4 

7 Q= J Kl.cos(w.t).dt + JAT2.^ 

8 

9 Thus , 

10 ti + T/4 

11 Q = [ (Kl/w) .sin(w. t) ] + K2.T/4 

12 ti 

13 

14 or 
15 

16 Q = (Kl/w) . [sin (w.ti + w.T/4) - sin(w.ti)] + L 

17 

18 or 

19 Q » (Kl/w) • [sin (w.ti + 7t/2) - sin (w.ti)] + L 

2 0 or 

21 Q = (Kl/w) . [cos (w.ti) - sin (w.ti)] + L (ii) 

2 2 where L is a constant = K2.T/4. 
23 

24 Using equation (ii) , the integration of V(t) from an 

25 arbitrary time tO to time tO+T/4 yields 
2 6 

27 Ql = (Kl/w) . [cos (w. to) - sin (w. to)] + L (iii) 

28 

29 Using equation (ii) , the integration of V(t) from 

30 time tO+T/4 to time tO+T/2 yields 
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1 Q2 = (Kl/w) . [cos (w.tO +W.T/4) - sin(w.tO + w.T/4)]+L 

2 or 

3 Q2 = (Kl/w) . [cos (w. tO + 7i/2) - sin(w.fcO + 7t/2)]+L 

4 or 

5 Q2 = (Kl/w) . [-sin (w. tO) - cos(w.tO)] + L . ..(iv) 

6 

7 Using equation (ii) , the integration of V(t) from 

8 time tO+T/2 to tO+3T/4 yields 
9 

10 Q3 = (Kl/w) . [cos (w.tO+w.T/2) - sin (w. tO+w. T/2 ) ] +L 

11 or 

12 Q3 = (Kl/w) . [cos (w. t0+7i) - sin (w. t0+7t) ] + L 

13 or 

14 Q3 = (Kl/w) . [-cos (w. tO) + sin (w. tO)] + L • (v) 

15 

16 Using equation (ii) , the integration of V(t) from 

17 time tO+3T^4 to time tO+T yields 
18 

19 Q4 = (Kl/w) . [cos (w.tO+w.3T/4) - sin (w. tO+w. 3T/4) ] +L 

2 0 or 

21 Q4 = (Kl/w) . [cos (w.t0 + 37i/2) - sin (w . 1 0 + 371/2 ) ] +L 

2 2 or 

23 Q4 = Kl/w) . [sin (w. tO) + cos(w.tO)} + L (vi) 

24 

25 Writing K = Kl/w and a = w.tO, then equations (iii) 

2 6 through (vi) yield for the four successive 

27 integrations of V(t) 
28 

29 - Ql = -K. since + K.coscc + L (vii) 

30 Q2 = -K.sina - K.cosa + L • (viii) 

31 Q3 = K.sina - K.cosa + L (ix) 



Q4 = K.sina + K.cosa + L 



(x) 



Integration control 

In order to control the sensor output integration, 
as just described, over four successive quarter 
periods of the drill string rotation, a train of n 
(with n any multiple of 4) equally spaced pulses per 
revolution must be generated. If one pulse P 0 of 
this pulse train is arbitrarily chosen at some time 
tO, the repeated pulses P n /4, Pn/2 and P 3n /4 define 
times tO+T/4, tO+T/2 and tO+3T/4 respectively where 
the period of rotation T = 2tt/w and w is the angular 
velocity of rotation. 

A suitable means for generating an appropriate 
control pulse train is described in US-A1- 
20020078745, which is hereby incorporated by 
reference. 

In an alternative form of integration control, the 
sensor output waveform itself can be used with 
appropriate circuitry for defining the integration 
quadrant periods. In particular, the relatively low 
noise magnetic fluxgate output is well suited to act 
as input to a phase -locked- loop arrangement. Fig. 2 
shows such an arrangement, successive output pulses 
defining the integration quadrants. 



Rotation angle 
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1 Equations - (vii) through (x) can be solved to yield 

2 angle a; there is a degree of redundancy in the 

3 . possible solutions but, for example, 
4 

5 Ql - Q2 = 2K.coscc 

6 and 

7 Q3 - Q2 = 2K. since 

8 or 

9 sina/cosa = (Q3 -Q2 ) / (Q1-Q2 ) (xi) 

10 

11 Since a= w.tO, the angle S(t0) between the axis 

12 (OX) and the direction of {VOXY} at time tO can be 

13 determined from equation (xi) , and the angle between 

14 (OX) and {VOXY} at any time tm measured from the 

15 arbitrary starting time tO is then 
16 

17 S(tm) = a + w.tm = S(t0) + 27C.tm/T . . . . (xii) 

18 

19 Magnitudes of vectors {VOXY} and {VOZ> 

20 

21 Equations (vii) through (x) can be solved to yield 

22 the constant L: 
23 

24 L = (Ql + Q2 + Q3 + Q4)/4 (xiii) 

25 

2 6 and the constant K can be determined from: 
27 

28 (K) 2 = [(Ql-L) 2 + (Q2-L) 2 ]/2 

29 = [(Q3-L) 2 +(Q4-L) 2 ]/2 . . . (xiv). 

3 0 

31 The magnitude of vector {VOZ} can be determined as 
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18 
19 



26 
27 
28 
29 
30 



1 

2 VOZ = (K2-c)/cos(m) = (4.L/T - c)/cos(m) (xv) 

3 provided that constant c is known. 
4 

5 The magnitude of vector {VOXY} can be determined as 
6 

7 VOXY = Kl/sin(m) = (K.w) / sin(m) (xvi) 

8 

9 Inclination angle 
10 

11 The inclination angle (INC) can be derived from the 

12 gravity vector {G} with the aid of a rotating 

13 accelerometer. 
14 

15 Referring to Fig. 3, where (INC) is the angle 

16 between the tool axis (OZ) and the gravity vector 
IV {G>, 



GOZ = G. cos (INC) (xvii) 

2 0 and 

21 GOXY = -Gsin(INC) (xviii) 

22 

2 3 The accelerometer output can be written as 
24 

2 5 VG(t) = GOZ. cos (m) + GOXY . sin (m) . cos (wt) 

+ CP. sin (m) + D . sin(m) (xix) 



where CP is a centripetal acceleration term and D is 
a sensor datum term. The centripetal acceleration 
term CP is zero for configuration. 2 and makes this 

31 the preferred configuration for mounting of the 

32 accelerometer. 
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'1 

2 Since CP is proportional to w 2 /r and is constant for 

3 constant w, then clearly VG(t) is of the form 
4 

5 VG(t) = Kl.cos(w.t) + K2 (w) 

6 . (or Kl.cos(w.t) + K2 for configuration 2) (xx) 

7 

8 where Kl and K2 (w) are constants at constant angular 

9 velocity w in the case of configuration 1 and always 

10 constant in the case of configuration 2. the 

11 constants Kl and K2 (w) can be determined from the 

12 accelerometer output integrations as described above 

13 together with the angle (Highside Angle HS = w.t) 

14 between the axis (OX) and the direction of {GOXY} . 
15 

16 Kl = GOXY. sin (m) (xxi) 

17 and 

18 K2 (w) = GOZ.cos(m) + D.sin(m) (xxii) 

19 with. 

20 C(w) = CP.sin(m) + D.sin(m) (xxiii) 

21 constant at constant angular velocity w (or for 

22 configuration 2 at all w) . 
23 

24 A calibration procedure can be carried out to 

2 5 determine the values of C (w) for angular velocity 

2 6 values w (constant in the case of configuration 2) 

27 by calculating values of K2 (w) with the rotation 

28 axis (OZ) horizontal when C (w) = K2 (w) . 
29 

3 0 Thus, for any drilling situation with known angular 

♦ 

31 velocity w, the vector components of the local 

32 gravity vector {g} can be determined as 



4 



1 

2 GOXY = Kl/sin (m) - ....... (xxiv) 

3 and 

4 GOZ = (K2 (w) - C (w) ) /cos (m) (xxv) 

5 



6 The inclination angle (INC) can then be determined 

7 from 
8 

9 sin (INC) /cos (INC) = -GOXY/GOZ (xxvi) 

10 

11 . Azimuth angle 
12 

13 When using a rotating fluxgate, the azimuth angle 

14 (AZ) can be determined from a consideration of the 

15 magnetic vector {b} . What follows is applicable to 

16 both configuration 1 and configuration 2. 
17 

18 With reference to Fig. 4, it can be shown that 
19 

2 0 BOZ = BV. cos (INC) 

21 + BN.cos (AZ) . sin (INC) (xxvii) 

22 

23 and 
24 

25 BOXY = (BN.cos (AZ) .cos (INC) -BV. sin (INC) ) .cos (HS-MS) 

26 + BN. sin (AZ) .sin (HS-MS) (xxviii) 

27 

2 8 or, with HS-MS = d a constant ,< 
29 

3 0 BOXY = (BN.cos (AZ) . cos (INC) -BV. sin (INC) ). cos (d) 

31 +BN.sin(AZ) .sin(d) (xxix) 

32' 
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13 

1 • With D the fluxgate datum, the fluxgate output can 

2 be written 
3 

4 VB(t) = BOZ. cos (m) + BOXY. sin (m) .cos(w.t) 



5 + D.sin (m) (xxx) 

6 or 

7 VB(t) = Kl.cos (w.t) + K2 (xxxi) 

8 where 

9 Kl = BOXY • sin (m) 

10 and 



11 K2 = BOZ.cos(m) + D.sin(m) 

12 = BOZ.cos(m) + C (xxxii) 

13 

14 are constants which can be determined from the 

15 fluxgate output integrations as described above 

16 together with the angle (Magnetic Steering Angle = 

17 MS = w.t) between the axis (OX) and the direction of 

18 {BOXY}. 
19 

2 0 A calibration procedure can be carried out to 

21 determine the value of the constant C by calculating 

22 the value of K2 while rotating about the direction 

23 of the axis (OZ) along which BOZ = 0 when K2 = C . 
24 

25 Thus, for any drilling situation the vector 

26 components of the local magnetic field {b} can be 



2 7 determined as 
28 

29 BOXY = Kl/sin (mj (xxxiii) 

3 0 and 

31 BOZ = (K2-C)/cos(m) (xxxiv) 

32 
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With reference to Fig. 5, the horizontal component 
{BN} of the local magnetic field vector {b} can be 
represented by horizontal components {Bl} and {B2 } 
where 

Bl = BOXY. cos (d) .cos (INC) 

+ BOZ . sin (INC) (xxxv) 

and . 

B2 = BOXY, sin (d) (xxxvi) 

The Azimuth Angle (AZ) can then be determined from 

sin(AZ) /cos (AZ) = -B2/B1 (xxxvii) 

Also, the horizontal component of the local magnet: 
field can be determined from • - 

BN = (Bl 2 + B2 2 ) 3/2 (xxxviii) 

and the vertical component of the local magnetic 
field can be determined from 

BV = BOZ. cos (INC) 

- BOXY. cos (d) .sin (INC) (xxxix) 

Earth's rotation vector 



Where it is not practicable to use a magnetic 
fluxgate, this may be replaced by a rate gyro as 
sensor. 
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1 With reference to Fig. 6, if the geographic latitude 

2 at the drilling location is (LAT) then the vertical 

3 component of the earth's Rotation Vector {RE} is 



4 

5 RV = -RE.sin(LAT) (xl) 

6 and the horizontal component is 

7 RN = RE.cos(LAT) (xli) 

8 



9 The magnitude of the cross -axis rate vector {ROXY} 

10 can be shown to be 
11 

12 ROXY = (RN. cos (GAZ) . cos (INC) -RV. sin (INC) ) . cos (d) 

13 + RN. sin (GAZ) sin (d) (xlii) 

14 

15 where (GAZ) is the gyro azimuth angle and 

16 d = HS - GS is constant . 
17 

18 Since RN # RV, d and INC are known and ROXY can be 

19 derived as discussed below, (GAZ) can be determined. 
20 

21 With the particular configuration where the rate 

22 gyro sensing axis is perpendicular to the drill * 

23 string rotation axis (OZ) , the -rate gyro output can 



24 be written 
25 

26 VG(t) = ROXY. cos (w.t) + D (xliii) 

27 

2 8 where D is the rate gyro datum, or 
29 

30 VG(t) = Kl.cos(w.t) + K2 (xliv) 

31 



16 

1 where the constant Kl = ROXY can be determined from 

2 the rate gyro output integrations as described above 

3 together with the Gyro Steering Angle GS = w.t 

4 between (OX) and the direction of {ROXY} . 

' 5 

6 The variation in the Rate Gyro Datum makes it 

7 difficult to achieve satisfactory datum calibration 

8 in all circumstances. It is unlikely that Gyro 

9 Azimuth measurements should be attempted at high 

10 inclination angles. The use of the rate gyro is 

11 most likely with near-vertical boreholes in 

12 locations where magnetic azimuth measurements are 

13 unreliable (such as close to rigs) and the Gyro . 

14 Azimuth GAZ is approximately equal to the angle d. 
15 

16 The present invention thus makes possible the 

17 measurement of a number of borehole-related 

18 parameters during rotation of a drillstring and 

19 using a reduced number of sensors. Modifications 

2 0 may be made to the foregoing embodiments . within the 

21 scope of the present invention. 




V(t) ■ Kl.cos(w.t) + K2 
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